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Abstract 

YBa2Cu3C>7 (YBCO) exhibits a large anisotropy between the a and b axes 
in the Cu02 planes because of the presence of CuO chains. In order to ac- 
count for such an anisotropy we develop a Ginzburg-Landau (GL) theory for 
an anisotropic d-wave superconductor in an external magnetic field, based 
on an anisotropic effective mass approximation within CuC>2 planes. The 
anisotropic parameter A = m x /m y} where m x (m y ) is the effective mass in 
the x (y) direction, is found to have significant physical consequences: In the 
bulk case, there exist both the s- and d-wave order parameters with the same 
transition temperature, as long as A / 1. The GL equations are also solved 
both analytically and numerically for the vortex state, and it is shown that 
both the s- and d-wave components show a two-fold symmetry, in contrast to 
the four-fold symmetry around the vortex, as expected for the purely d-wave 
vortex. With the deviation of A from unity, the opposite winding between 
the s- and d-wave components observed in the purely d-wave case is gradu- 
ally taken over by the same winding number. The vortex lattice is found to 
have oblique structure in a wide temperature range with the precise shape 
depending on the anisotropy. 

PACS numbers: 74.20.De, 74.60.-w, 74.60.Ec, 74.72.-h 
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I. INTRODUCTION 



Recently, the order parameter symmetry has become the central issue in high-T c super- 
conductivity Many experiments which probe directly the phase of the pairing state have 
provided strong evidence for a sign change of the order parameter |]-§]], consistent with a 
predominantly <i-wave pairing symmetry. At the same time, there are several measurements 
which can not be explained within the simple d x 2_ y 2 state. For example, finite tunneling 
current along the c axis of a copper oxide clearly shows an s-wave character |J, because 
there should be no Josephson current between a Cu0 2 plane with a gap of d x 2_ y 2 symmetry 
and a conventional s-wave superconductor. 

It is well known that YBCO is not in the purely tetragonal phase due to the existence 
of chains. Indeed, YBCO exhibits a large anisotropy between the a and b directions in the 
measurements of the penetration depth || and the vortex structure by scanning tunneling 
microscopy (STM) [[7J. It was argued recently || that these apparently conflicting exper- 
imental results in this material may be explained by assuming that there exist two order 
parameters, with different symmetry but the same transition temperature. Namely, the 
main gap with a <i-wave symmetry would result from the Cu02 planes, and a smaller s-wave 
component would be due to the CuO chains. 

In this work, we will consider a simple model for an anisotropic d-w&ve superconductor, 
based on the anisotropic effective mass approximation within a single Cu02 plane. In this 
model, the a-b anisotropy of YBCO is taken into account by a single parameter, namely 
the electron mass anisotropy, A = m x /m y , which can be fit to the measured penetration 
depth anisotropy ||. Then the Ginzburg-Landau (GL) theory for such an anisotropic d- 
wave superconductor will be studied. Following the procedure described in Ref . || , we first 
derive microscopically the GL equations for this anisotropic system, and then consider the 
possible solutions of these GL equations for both bulk and vortex states. We will show 
that the anisotropic parameter A has significant physical consequences: in the bulk, the 
stable solution from our theory is the mixed s + d state, and both the s- and d-wave order 
parameters have the same transition temperature. This s + d state is just what we want to 
explain the tunneling data and other apparently conflicting results observed in YBCO. The 
GL equations are also solved both analytically and numerically for the vortex structures. 
We find that the anisotropic <i-wave vortex is very different from the purely (i-wave case. 
Namely, both the s- and (i-wave components show a two-fold symmetry, in contrast to the 
four-fold symmetry around the vortex as expected for the purely ci-wave vortex. Specifically, 
the (i-wave order parameter exhibits an elliptic shape and the s-wave component shows a 
shape of butterfly. With the deviation of A from unity, the opposite winding between the s- 
and (i-wave components obtained in the purely <i-wave case is gradually taken over by the 
same winding number. The vortex lattice is found to be in oblique in a wide temperature 
range with the precise shape depending on the anisotropic parameter A. Here we wish to 
point out that the preliminary results of A = 1 case for the structures of a single vortex 
and vortex lattice were reported recently in a conference on superconductivity [ID]. Here we 
concentrate on the results for the A ^ 1 case. 

In Sec. II, starting from Gorkov's theory of weakly coupled superconductors JO], the 
GL equations for the anisotropic ci- wave superconductor are derived. In Sec. Ill, we discuss 
the possible solutions of the GL equations for a uniform or bulk system. In Sec. IV, we 



2 



study analytically the qualitative features of a single vortex using the GL equations. In 
Sec.V, we present the numerical result for single vortex structure. The numerical calculation 
for a vortex lattice structure is performed in Sec. VI, and Sec. VII includes conclusion and 
discussions. 



II. GL EQUATIONS 



In this section we shall derive the GL equations for the anisotropic d-w&ve superconduc- 
tor, following closely the procedure we used for a purely d-wave superconductor ||. Here 
only the main steps are presented and the difference between the present work and the 
previous one will be ephasized. Our starting point is the gap equation 



A*(x, x') = V(x - x')T ^(x, x', u n ) 



(2.1) 



which allows for more general than conventional s-wave pairing. V(x — x') is the effec- 
tive two-body interaction of the weak-coupling theory. Using Gorkov description of a 
superconductor in the magnetic field: 

%<jj n - 77— (— «V + eA) 2 + /i 
2m 



—ioo n - J— (zV + eA) 
2m 



G(x, x', uj n ) + J dx"A(x, x")F+(x", x', uj n ) = 5(x - x'), (2.2) 
F+(x, x', w n ) + / dx"A*(x, x")G(x", x', u n ) = 0, (2.3) 



where // is the Fermi energy and A is the vector potential. The normal-state Green's function 
at zero magnetic field can be written in the form 



Gq(x, LO Ti 



(27T) 



c/ke 



ikx 



where 



+ 



2m x 2m y 



ILOr, 



(2.4) 



(2.5) 



is the single particle energy measured from the Fermi energy fi. We note that the case 
with m x = m y in the above equation corresponds to the isotropic d-wave superconductor. 
Keeping up to the third order in A, defining the central-mass coordinates R = (x + x')/2 
and the relative coordinates r = x — x', and making Fourier transform with respect to the 
relative coordinates, we obtain 01 



A(R, k) 



dW 1 

nk'-k)TE^rA(R,k') 



(2tt)2 
dk' 
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A(R, k') 
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2 \2 



A(R,k')| 2 A(R,k'), 



(2.6) 
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where we have introduced the operator 



n = iVj 



2eA 



R- 



(2.7) 



In order to obtain the generic Ginzburg-Landau equations, which govern the spatial vari- 
ation of the order parameters, for an anisotropic ci-wave superconductor, we need to specify 
the form of the interaction. Here we use the following ansatz for the effective interaction 
responsible for the spin-singlet pairing: 



y( k - k') = -V a + V d (k 2 x - k 2 y )(k' 2 - ky), 



(2i 



where k = k/|k| is the unit vector in the direction of k. By taking both Vd and V s positive, 
the —Va corresponds to attractive interaction responsible for <i-wave pairing, and V s can 
be regarded as an effective on-site repulsive interaction. The general expression of order 
parameter that follows Eq.(2.8) is 



A(R,k) = A s (R) + A d (R)(^_^). 



(2.9) 



Substituting Eqs.(2.8) and (2.9) into Eq.(2.6), and Comparing both sides of the gap equation 
for /^-independent terms and terms proportional to (k 2 — k 2 ), we obtain the GL equations 
in a form suitable for finding the GL free energy functional: 



+ 



A - 1 



A + l 
A + 2V 7 ! 



A d + 2«7 d /i< 



n 2 



+ 



^ 2m x 2m y / 



1 IT; 2\/A + 1 



A II 2 



(1 + V\) 2 2m, (1 + V\) 2 2m y 



A ^ + f+>' 2A * + ^ A ^ 



+ l (l + v/A) 2 



AJAJ 2 + 



M/A - 1 

(TTTaY 3 



0, 



(2.10) 



1 - Id hi 



(2e<uj Q \ 2(1 + A) 



V ttT J (1 + v/A)2j A + l 



A„- 



+2«7d/i 



1 - >/A + 3A + \V\ II 2 1 + 3^ - A + Av^ n 2 



;i + Vxy 



2m 7 



;i + V\y 



2m,, 



A, 



+ 



a + 2V\ - 1 n 2 2V\ + 1 - a n 2 



+27d« 



l + v^)2 2m 

Vx-i 



+ 



v/A + l 
A + l 



AfA s * + 2- 



T + v^)2 2m ?/ 



A, 



i + vxy 



A 2 A* + +2 



-v/A) 2 
Av^- 1 

(TT7a> 



a, + ^^La:a 2 



s d 



A,|A,| 2 + 



(1 + v/A) 3 
T + A)(l + A + v / A) 
(1 + v/A) 4 



A d | 2 A c 



0, (2.11) 



where a = 7((3)/8(irT) 2 , 7 is the Euler constant, 7 d = N(0)V d /2 is the interaction strength 
in the purely <i-wave channel when A = 1, and 
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(i + Va) 2 ^ , (i + Vx) 2 v s 

as ~ 1 + A l i+ 1 + A V d 



(2.12) 



It is easy to show that if setting A = 1, Eqs.(2.10) and (2.11) return back to the results we 
obtained for a purely <i-wave superconductor ||. The corresponding GL free energy is 



1 - Id In 
+27 d a/x 



(2e 1 u D \ 2(1 + A) 



V ttT J (1 + v^) 2 
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+ h.c. 
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(1 + v/A) 2 ' ' 
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;i + Vxy 



|A d | 2 (A:A, + A s A*) 



(2.13) 



It is interesting to note in the above equation that except for the mixed gradient terms, 
which are induced by the magnetic field, there exist the new terms, such as 



oc(A-l)(A:A d + A s A^) 



and 



oc 



(Va-i)( 



AJ 2 + 1 + v/X+A 



A d \ 2 Wa:a, + a s a: 



These new terms come completely from the mass anisotropy. For isotropic systems with 
A = 1, these terms vanish. In the following section, we will discuss the physical consequences 
of these new terms. 



III. BULK SOLUTIONS 



In this section we study the solutions of the GL equations of an anisotropic c?-wave 
superconductor for a bulk or uniform system. In this case the gradient terms in Eqs.(2.10) 
and (2.11) are equal to zero. First let us examine the T c formula. For T — > T c , the coefficients 
of the linear terms in the GL equations determine the transition temperature: 




(i + Vxf 

2(1 + A) 




(3.1) 
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where T c0 is the transition temperature for a purely o?-wave superconductor corresponding 
to A = 1, which is given by 



In 



2e<uj D 1 



For a small anisotropy (A 



7tT c0 Id 

1), Eq.(3.1) reduces to 



(3.2) 



In 



T, 



1 



cO 



■(A-l) 2 >0. 



(3.3) 



This result implies that T c increases as the system deviates from the isotropy. 

In the following discussion, for convenience, we put the GL free energy into a dimen- 
sionless form, which can be done by scaling the energy by 4 7 d(l — T/T c ) 2 /3a, lengths by 



lfm/2m x (l-T/T c ), and setting A s = ^ S A~(A = 1), A d = i/i d Af(\ = 1), where 
A^A = 1) = 



F = 



4(1 - T/T c )/3a, and A = (2n£ x /$ )A: 

2(1 + A) r 1 / A-l > 

(1 + V\)n a sld (l-T/T c ) V A + 1, 
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1 - v / A + 3A + Av^ 
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7d (l - T/T c ) 
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Ajn^Vs 



+ 2 



(1 + v/A) 3 
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(1 + v/A) 2 



n^dl 2 + a 



1 + 3V\- A + A^A 



(l + V\f 

U y ip s U ip d + h.c. 



+ 



+ 



4(A + 1)(A + VX + 1) 4 + , 



3 (1 + V 7 !) 4 

4 A + 1 



1 + v/A) 2 "^ S " r 
4, |4 16 A + l 
3 IW ' 3 (1 + v/A) 2 



3(l + v / A) 2 
8 AVA- 1 



(3.4) 



In the above expression, the magnetic field energy has been included explicitly with k being 
the GL parameter. From 5F/5ip* and SF/Sip d we can obtain the GL equations for ip s and ijj d . 
We now discuss the bulk solutions. Assuming ip s = \ip s \e l0s , ip d = \ip d \e t9d , and 6 = 9 S — 9 d , 
we can determine the value of 6 through dF/d6 = and d 2 F/d6 2 > 0. We find that the 
stable solution is only possible for = 0, which means that the bulk system is in the mixed 
s + d state with a real combination. For T — > T c , ip s and ip d are given by 



*- = h (xtt) +* 



(3.5) 



(3.6) 
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with 



D(X) 



4_(l + yX) 2 [(A 
3" 



1)(A + VA + 1 




A-r 

A + l. 



(3.7) 



It is easy to see that -D(A) — > 1 as A — > 1. Namely for an isotropic system with A — > 1, 
ipd — > 1, and — > 0, which implies that the purely d-wave state is only possible for the 
isotropic system. Whenever the system has an anisotropy, the mixed s + d state with if>^ < 1 
and ip s 7^ is generated. It is clear from Eq.(3.5) that these two order parameters have the 
same T c , which is given by Eq.(3.1). As argued by many authors ||, such a mixed s + d 
state in the bulk is just what we need to explain the tunneling data and other apparently 
conflicting results observed in YBCO. 



IV. VORTEX SOLUTIONS 

In this section, we determine the single vortex solutions for an anisotropic d-w&ve super- 
conductor using the GL equations derived in Sec. II. Previously, we have studied the purely 
ci-wave vortex structure |9j and found that near the vortex core, there coexist the s-wave 
and d-w&ve components with the opposite winding numbers. Far away from the vortex core, 
the induced s-wave component shows strong four-fold anisotropy and decays as r~ 2 . We 
expect that the mass-anisotropy will affect the vortex structure. For simplification, here we 
study the case when e = A — 1 is a small parameter. In this case the GL equations become 

a e 8 8 4 

7d(l - T/T c ) 27(1 - T/T c ) 3 3 3 

+\e{\M 2 ^d + V>M + \eMl>d + 2II 2 ^ + (IT 2 - IL> d = 0, (4.1) 



-^4 



;ipd - 



2 7d (l - T/T c 



4a s7 (l - T/T c ) 

+\zVI>.N. + \<^*s + \M 2 ^s) + n 2 ^ + (n 2 - n 2 ,)^ = 



^ + \H 2 *Pd + *\i>s\ 2 ipd + \^* d 



(4.2) 



In terms of the cylindrical coordinates, R = (r, 9), we expect that the c?-wave component 
has the form ip^ = e ld in the region of 1 r London penetration depth. Also note that, 
in this region, the magnetic field effect can be neglected. Then the leading terms in the 
equation (4.1) for ip s are 



a s , _1 
7d (l -T/T c f s 2 C 



7d(l - T/T c 



+U™r s - & - dl) e» = 0. 



(4.3) 



This equation suggests the following solution: 
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r 

where 

3e(3a s + 4) 



i) s = ae ie + \ (be- w - ce 3ie ), (4.4) 



(3a s + 8) 2 - 16 ' 



(4.5) 



b=* f 3 ^ + 2 ° , (4.6) 

2 (3a s + 8) 2 -16' V ; 

3 9a, + 28 

C 2 (3a, + 8) 2 -16' 1 j 

with a s = a s /jd(l — T/T c ) and e = |e[l/7d(l — T/T c ) — 1]. For T — > T c , ^ s takes the simple 
expression: 

_e_ e * + 7d(l ~T/T C ) 1 L _ lQ ^ m 
2a q 2a* r 2 



^ s = —e ty + Iax n ' (e- M - 3e My ) . (4 



It is very important to note that the first term in the above equation is independent of 
both temperature and the distance from the vortex core. This term comes, in fact, from 
the contribution of the bulk. Comparing this term with the bulk solution given in (3.4), we 
immediately find that they are identical for small e parameter. Our solution (4.4) implies 
that far away from the vortex core and the temperature approaches to T c , the bulk term 
of the s-wave component, with the same winding with respect to <i-wave order parameter, 
becomes dominant. The magnitude of ip s is 

|^ s |2 = a 2 + J_ ( fe 2 + c 2 } + 2a (6 _ c) cos20 _ Zbc cosAd 

This result clearly shows a two-fold symmetry due to the existence of the cos 2$ term. If 
setting e = 0, the cos 2(9 term vinishes, \ip s \ 2 recovers the four-fold symmetry, and our result 
returns back to that for a purely ci-wave superconductor, as given in Ref. 0. 
Near the vortex core, to the leading order, our GL equations become 

-ipd~ V 2 V rf = 0, (4.10) 



a. 



i> s - € _ i> d - 2V 2 ^ S - (dl - d 2 y ) fa = 0. (4.11) 



ld (l-T/T c y b 2 7d (l -T/T ( 
From (4.10) we have 

1 



fa = c (r - -r 3 j e l \ (4.12) 
where c is a constant. Putting the above equation into (4.11), we obtain 
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(4.13) 



Thus, the leading order terms of the order parameters near the vortex core are 



Core 
c r 

2a, 



ee w - 7d (l-T/T c 



-w 



(4.14) 
(4.15) 



For isotropic systems with e = 0, the above results reduce to those for a purely <i-wave 
superconductor. Namely, the s-wave component, with the opposite winding relative to the 
c?-wave order parameter, is induced near the vortex core ||. However, the anisotropy alters 
such a picture: As T — > T c , the opposite winding of the s-wave component is gradually 



taken over by the same winding term. Also e td and e" 
a two-fold symmetry around the vortex core. 



terms in Eq.(4.15) combine to give 



V. NUMERICAL RESULTS OF SINGLE VORTEX 



In last section, we have discussed analytically the asymptotic behavior of the single vortex 
for an anisotropic <i-wave superconductor. But the precise shape of the vortex structure is 
still not clear and it has to rely on the numerical calculation. Here, we perform a numerical 



study of the discretized GL free energy (3.4) using numerical relaxation approach ||12| , |T3 



In order to minimize the GL free energy functional in the presence of magnetic field, we use 
the constraint of fixing the average magnetic induction B by specifying the total flux $ in 
the unit cell, and impose the so-called "magnetic periodic boundary conditions" |13| , |i~4| . 

To perform the numerical relaxation calculation, we need to discretize the GL free energy 
(3.4) first. With the use of the forward difference approximation for the derivatives and 
taking into account the gauge invariance, we can write the free energy (3.4) in the discrete 
form: 



F — Fq + Fkin + Ffield, 



(5.1) 
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3 (i + Vxy 

8y/X- 1 



3^+1 
8 AVA- 1 
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3 (1 + VA) 



Mh j)\ 2 («, j)Mh j) + ^* d (h j)] 

\i>d(i,j)\ 2 [*P* s (hj)ipd(i,j) +ip d (ij)ip s (i,j)} 



(5.2) 
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ia x A x (i,j)\2 j fl 2 



+ X\,l; s (iJ + l)-i; s (i,j)e ia ^^\ 2 /a 2 y 



(i + V\) 3 



+A 



\MiJ + l)-MiJ)^ Av ™\V<i 



+ 



,ia x A x (i,j) 



(1 + VA) 3 
X + 2VX-1 
(1 + v/A) 2 

^Txy^ (V>.(<. i + 1) - Mi, 3Y ay 

x j + 1) - V£(«, j)^^) /a 2 , + h.c. 



i0x^a:(ij)l /„2 



(5.3) 



r 2 

"field 



K 2 f 1 2 

TTTf E + U) - - [^(i, j + 1) - A x (i,j)]/a y , (5.4) 

where (N y ) is the number of lattice points in the x (y) direction. On each lattice point 
the order parameters ip s and ipd have the values tp s (hj) and if)d{hj)> an d each point 

is associated with horizontal and vertical bonds. a x and a y are the lattice constants and 

A x (i,j) and A y (i,j) are the vector potential components on bonds — > (i + 1, j)] and 
— > + 1)], respectively. It is easy to show that in these lattice notations, the above 

expressions are invariant with respect to gauge transformation: 

A x (i,j) -> A x (i,j) + [x(i + - x(hj)]/a x , 

where x{hj) is the arbitrary phase of the order parameters at site Accordingly, 
the free energy and other physical quantities are also gauge invariant. To obtain simple 
boundary conditions we can choose a gauge such that A x is independent of x. In this case, 



our boundary conditions are |14]|L3" 



A x (0) = A x (L y ), (5.5) 
A y (L x ,y)-A y (0,y) = $/L y , (5.6) 
A y (x,Ly)=A y (x,0), (5.7) 

ips,d(x,L y ) = VvMK* /2 , (5.8) 
^ d (L XJ y)=^ d (0,y)e l ^ L y J (5.9) 

where L x = N x a x and L y = N x a y . In order to study the single vortex structure, we can 
choose one quantum of flux (i.e., $ = 2n) in a square unit cell with N x N lattice points. 
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With Eqs.(5.1)-(5.4) and the above boundary conditions, we can now realize the relaxation 
procedure. Choosing if> s , if)*, ip d , ifj* d , A x , and A y as independent variables, we can write 
down the relaxation iteration equations: 



^ n+1) ( 




= 4 n) ( 


hi) 


^ n+1) ( 


ij) 


= #( 


hi) 


A(n+1) 


(iJ) 


— A^> 




Ain+1) 






(i,3 



dF 



dip*(i,j) 
dF 



'dip* d (i,j) 
dF 



9A x (i,j) 
dF 



(n) 

j 

in) 

■> 

(n) 
(«) 



dAy(i,j) 



(5.10) 
(5.11) 
(5.12) 
(5.13) 



where e's are all positive numbers to be adjusted to optimize the convergence rate and n 
is an integer labeling the generations of iteration. It has been shown mathematically that 
F will monotonically decrease to its optimum state as n increases as long as we choose a 
proper initial state |I2| . 

In our numerical calculation, the parameters chosen are N x = N y = 101, a x = a y = 
0.2^, k = 2, T = 0.5T C , and V s = 0. With these parameters, the external magnetic 
field corresponds approximately to the thermal critical field H c . The use of the different 
parameters does not alter the qualitative physics. Let us first show the results for an 
isotropic <i-wave vortex structure with A = 1. Figs.l and 2 are typical surface plots for 
the distribution of the <i-wave order parameter and local magnetic field around the vortex, 
respectively, which look like the conventional s-wave vortex. But if looking at them closely, 
we find the difference from the conventional s-wave vortex. Fig. 3 is the contour plot of 
the o?-wave order parameter. It is clear that \if>d\ exhibits a four- fold symmetry. The local 
magnetic field also shows a similar four- fold anisotropy (not shown in the contour plot). 

The most interesting feature of a single vortex is that a small s-wave component is 
induced around the core, as shown in Fig.4 (A) (surface plot) and (B) (contour plot). One 
can clearly see that the distribution of \if) s \ exhibits the profile in the shape of a four- leafed 
clover, which is in agreement with our analytical result |J. We believe that the presence of 
this four-fold symmetric s-wave component is the reason to cause four-fold symmetry of the 
ci-wave order parameter \if) d \ (see Fig. 3) and the local magnetic field h around the vortex. 

We now discuss the anisotropic case with A / 1. With the deviation of A from unity, 
we find that both and h begin to show a two-fold symmetry, in contrast to the four- 
fold symmetry as expected for a purely <i-wave vortex. Specifically, both of them exhibit an 
elliptic shape. Fig.5 is the contour plot for \if>d\ with A = 2, which corresponds approximately 
to the measured panetration depth anisotropy . 

The s-wave component is much more sensitive to the anisotropic parameter A. Fig. 6 
(surface plot) and Fig. 7 (contour plot) show how changes with the increase of A. It 
is apparent that \if) s \ exhibits the two- fold symmetry, and its shape changes from the four- 
leafed clover in the purely (i-wave case (A = 1) to a butterfly as A increases. These results 
agree well with our analytical calculation [see Eqs.(4.8) and (4.15)]. 

Recently, the vortex structure of YBCO has been directly observed using STM imaging 
technique 0. An elongated shape of the vortex was realized. The ratio of the axes in 
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the apparent elliptic shape is about 1.5. Furthermore, this elongation was found to be 
independent of the scanning direction of the STM tip. We believe that this elongation 
directly reflects the a-b anisotropy. This observed vortex shape can be qualitatively account 
for by the present GL theory for an anisotropic (i-wave superconductor (see Fig.3). 



VI. NUMERICAL RESULTS OF VORTEX LATTICE 

Recent observation of an oblique vortex lattice structure in YBCO has been reported 
with the angle between the primitive axes (3 ~ 73° by the neutron scattering 0] and 
(3 ~ 77° by STM measurements |7]. The rich and complicated structure of a single vortex 
in an anisotropic <i-wave superconductor obtained in the previous sections will be expected 
to form a different vortex lattice than the conventional s-wave superconductor, and may 
provide an explanation to the oblique vortex lattice structure observed in YBCO. To check 
this, the vortex lattice structure is going to be studied using numerical relaxation method. 
The vortex lattice structure is still described by the discrete GL free energy functional given 
in Eqs.(5.1)-(5.4). We chose a rectangular unit cell with two vortices [13]. The periodic 



boundary conditions, very similar to Eqs.(5.5)-(5.9) except for $ = Air in the present case, 
are used in our calculations. The ratio of a y /a x controls the shape of the vortex lattice 
structure. For example, a y /a x = 1 corresponds to the square, while a y /a x = \/3 corresponds 
to triangular lattice. We have calculated the dependence of the free energy on the ratio of 
a y /a x using the same set of parameters as for the single vortex. 

Let us first present the results for the isotropic systems. Fig. 8 displays the free energy 
as a function of a y /a x for the isotropic (i-wave superconductor (A = 1). It is evident from 
the figure that the minimum of the free energy is at the position with a y /a x ~ 1.3, signaling 
that an oblique vortex lattice with the angle (3 ~ 75° between the primitive axes is stable. 
In this case, the superconducting state in the bulk or uniform system is purely <i-wave. The 
s-wave component near the vortex core is induced due completely to the mixed gradient 
terms in the GL free energy. To correlate between the single vortex structure and the vortex 
lattice, we have performed the calculations for different temperatures. We find that in a 
wide temperature range below T c , the presence of a sizable induced s-wave component causes 
four-fold symmetry of the (i-wave order parameter and the local magnetic field around the 
vortex. Such an anisotropic single vortex tend to form an oblique lattice structure. Fig. 9 
shows the oblique vortex lattice formed at T/T c = 0.5 by the (i-wave order parameter (A) 
and the s-wave component (B) for the isotropic case (A = 1). The local magnetic field 
distribution is very similar to the (i-wave order parameter (not shown in the figure). 

However, when T — > T c , the induced s-wave component is strongly suppressed, which can 
be also seen from our analytical results given in Eqs.(4.8) and (4.13). In this case, the induced 
s-wave component is too small to affect the distribution of the (i-wave order parameter and 
the local magnetic field. Consequently, both ipd and h have isotropic distribution around the 
vortex core, similar to the conventional s-wave vortex. These isotropic isolated vortices prefer 
to have a triangular vortex lattice, identical to the vortex lattice in an s-wave superconductor 



iql , as expected. This result is confirmed by our numerical calculation for the free energy 
that the minimum of the free energy moves to a y /a x = i.e., the triangular lattice is 
stablized. 

We now turn to the anisotropic case. Recent STM measurements in YBCO reveaved 
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an oblique vortex lattice with (3 ~ 77° [[/]]. Moreover, the elongated vortex cores with the 
ratio of principle axes about 1.5 were also found by this technique. As noted by the authors 
of Ref. 0, if this elongation reflects the intrinsic a-b anisotropy in a conventional s-wave 
superconductor, such an anisotropy would lead to a distorted vortex lattice with an angle 
(3 = 82° inconsistent with the observed value. Thus it seems that the a-b anisotropy alone 
can not explain the observed vortex lattice structure and additional effects, such as the order 
parameter symmetry, must be involved in order to account for the experimental data. We 
should note also that although the vortex lattice in a purely (i-wave superconductor has an 
angle |HJ[nj very close to the observed value, the four-fold symmetric vortex cores obtained 
in this case are inconsistent with the observed elliptic vortex cores. In this regard, it is inter- 
esting to study the vortex lattice structure in an anisotropic d-wave superconductor which 
contains both the a-b anisotropy as well as the (i-wave order parameter symmetry. Fig. 10 
shows the free energy as a function of a y /a x for A = 2, which corresponds approximately to 
the experimental data on the penetration depth |3j and the coherence length 0. It is clear 
that the minimum of the free energy locates at a y /a x ~ 1.3, almost the same position as for 
the isotropic (i-wave superconductor. 

We should mention that in an anisotropic (i-wave superconductor, the the a-b anisotropy 
and the order parameter symmetry play very different roles in determining the vortex struc- 
tures. The vortex tends to have a two-fold symmetry due to the a-b anisotropy, while it 
tends to have a four-fold symmetry due to the d-wave pairing state. For small A, the (i-wave 
order parameter symmetry is important, and the anisotropy becomes dominant for large 
anisotropy, as long as the effect on the vortex structure is concerned. We would like to 
point out here that the similar oblique vortex lattice obtained in both isotropic (A = 1) 
and anisotropic (A = 2) (i-wave superconductors is only a coincidence. In general, the real 
single vortex and the vortex lattice structures are determined by the competition between 
the anisotropy and the (i-wave order parameter symmetry. 

Fig. 11 shows the vortex lattice formed by the (i-wave (A) and s-wave (B) order parameters 
for an anisotropic (i-wave superconductor with A = 2. The local magnetic field distribution 
(not shown in the figure) is very similar to the (i-wave component. It is seen that the vortex 
lattice is oblique with j3 ~ 75° and the vortex cores are elliptic. These results are in perfect 
agreement with the STM measurements on the vortex lattice structure in YBCO f/J. 

VII. CONCLUSIONS 

We have established microscopically a GL theory for an anisotropic superconductor with 
d x 2_ y 2-wave pairing symmetry within the anisotropic effective mass approximation. The 
experimental basis of our model is recent measurements which revealed a large anisotropy of 
the penetration depth || and the coherence length between a and b directions in YBCO 
due to the existence of the CuO chains. In our model, a single anisotropic parameter A 
is introduced which gives a measure of the difference in the two effective masses m x and 
m y and this parameter can be fitted to the measured penetration depth and the coherence 
length data. The GL equations obtained in the present work should be useful to study the 
various properties of YBCO. 

We have considered the solution of the GL equations for a uniform or bulk system and 
found that the stable solution is the mixed s + d state, and both the s- and (i-wave order 
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parameters have the same transition temperature. Such an s + d state is just what we need 
to explain the tunneling data || and the other apparently conflicting experimental data 
observed in YBCO. 

We also solved the GL equations both analytically and numerically for the vortex struc- 
tures. For the single vortex we find that the anisotropic (i-wave vortex is very different from 
the purely (i-wave case. Namely, both the s- and (i-wave components show a two-fold sym- 
metry, in contrast to the four-fold symmetry around the vortex as expected for the purely 
(i-wave vortex. Specifically, the (i-wave order parameter exhibits an elliptic shape and the 
s-wave component shows a shape of butterfly. With the deviation of A from unity, the op- 
posite winding between the s- and d-wave components observed in the purely (i-wave case 
is gradually taken over by the same winding number. 

The vortex lattice is in general oblique for both the purely <i-wave and the anisotropic 
c?-wave superconductors. Although the angle between the primitive vectors of the vortex 
lattice in a purely c?-wave superconductor is comparable to the observed value, the shape 
of the vortex cores are very different from the experiments. On the other hand, for sn 
anisotropic (i-wave superconductor, the shape of the vortex lattice is determined by the 
competition between the anisotropy and the (i-wave order parameter symmetry. By using 
the anisotropic parameter obtained from the experimental data on the penetration depth and 
the coherence length, we were able to find a vortex lattice which agree well with experiments 
not only in angle between the primitive axes but also the elliptic shape of the vortex cores. 
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FIGURES 



Fig.l Surface plot for the distribution of |^| around a single vortex in the isotropic 
ci-wave superconductor. The parameters used are given in the text. 

Fig. 2 Surface plot for the distribution of the local magnetic field h around a single vortex 
in the isotropic rf-wave superconductor. 

Fig. 3 Contour plot for the distribution of around a single vortex in the isotropic 
d-wave superconductor. 

Fig. 4 Distribution of |^ s | around a single vortex in the isotropic rf-wave superconductor. 
(A) Surface plot, and (B) contour plot. 

Fig. 5 Contour plot for the distribution of \ipd\ around a single vortex in the anisotropic 
d-wave superconductor with A = 2. 

Fig. 6 Surface plot for the distribution of \ip s \ around a single vortex in the anisotropic 
d-wave superconductor with different anisotropic parameters: (A) A = 1.05, (B) A = 1.2, 
and (C) A = 2. 

Fig. 7 Contour plot for the distribution of \ip s \ around a single vortex in the anisotropic 
ci-wave superconductor with different anisotropic parameters: (A) A = 1.05, (B) A = 1.2, 
and (C) A = 2. 

Fig. 8 Free energy as a function of the ratio of a y /a x for vortex lattice in an isotropic 
d-wave superconductor. 

Fig. 9 Contour plots of the d- (A) and s-wave (B) order parameters for the stable vortex 
lattice structure in an isotropic rf-wave superconductor with a y /a x = 1.3, corresponding to 
an oblique lattice with (3 = 75°. 

Fig. 10 Free energy as a function of the ratio of a y /a x for vortex lattice in an anisotropic 
o?-wave superconductor with A = 2. 

Fig. 11 Contour plots of the d- (A) and s-wave (B) order parameters for the stable vortex 
lattice structure in an isotropic (i-wave superconductor (A = 2) with a y /a x = 1.3, corre- 
sponding to an oblique lattice with (3 = 75°. 
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